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MODULATIONAL INSTABILITY IN THE WHITHAM 
EQUATION WITH SUREACE TENSION AND VORTICITY 


VERA MIKYOUNG HUR AND MATHEW A. JOHNSON 


Abstract. We study modulational stability and instability in the Whitham 
equation, combining the dispersion relation of water waves and a nonlinearity 
of the shallow water equations, and modified to permit the effects of surface 
tension and constant vorticity. When the surface tension coefficient is large, we 
show that a periodic traveling wave of sufficiently small amplitude is unstable 
to long wavelength perturbations if the wave number is greater than a critical 
value, and stable otherwise, similarly to the Benjamin-Feir instability of gravity 
waves. In the case of weak surface tension, we find intervals of stable and 
unstable wave numbers, whose boundaries are associated with the extremum 
of the group velocity, the resonance between the first and second harmonics, the 
resonance between long and short waves, and a resonance between dispersion 
and the nonlinearity. For each constant vorticity, we show that a periodic 
traveling wave of sufficiently small amplitude is unstable if the wave number 
is greater than a critical value, and stable otherwise. Moreover it can be made 
stable for a sufficiently large vorticity. The results agree with those based upon 
numerical computations or formal multiple-scale expansions to the physical 
problem. 


1. Introduction 

The Korteweg-de Vries (KdV) equation 

( 1 . 1 ) Ut + ^/^(l+^cf^ljua;+UUa;=0 

approximates the water wave problem in a small amplitude and long wavelength 
regime, and furthermore, it satisfactorily explains solitary waves and cnoidal wave 
trains of the physical problem. Here, t G IR is proportional to elapsed time and 
x G M is the spatial variable in the predominant direction of wave propagation; 
u = u(x, t) is real-valued, related to the surface displacement from the undisturbed 
fluid depth d, and g is the constant due to gravitational acceleration. Throughout 
we express partial differentiation either by a subscript or using the symbol d. When 
waves are long compared to the fluid depth so that fed ^ 1, fc the wave number, 
one may expand the phase velocitjQfor the water wave problem in the irrotational 
setting and write that 

(1.2) cww(k) := ^ ^ + 0(fc^d"‘). 
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*The dispersion relation of water waves is well-known, whose derivation dates back to the work 


of Airy 


Air45 


in 1845! 
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Therefore the KdV equation may be regarded as to approximate up to second order 
the dispersion relation of the water wave problem in the long wavelength regime. 
Here and elsewhere, ± mean left and right moving waves. (One may take, without 
loss of generality, the + sign in (1.1). As a matter of fact, x i-A —x and u i-A —u 
turn one equation to the other.) 

However the KdV equation misses peaking into sharp crests. Furthermore waves 
in shallow water at times break into bores whereas the KdV equation prevents 
singularity formation from solutions. This is not surprising since the phase velocity 
associated with the linear part of the KdV equation poorly approximate^ cww 
when kd becomes large. One may incidentally argue that peaking and breaking are 
high frequency phenomena, for which short wave components become important 
and the long wavelength assumption is no longer adequate. 

Whitham therefore emphasized in [Whi74| that “it is intriguing to know what 
kind of simpler mathematical equation (than the physical problem) could include” 
breaking and peaking, and he put forward 

(1.3) ut + Mux + uux = 0, 

where is a Fourier multiplier, defined via its symbol as 


(1.4) 


Mf{k) = m{k)f{k), 


and m = cww (see (1.2|). It combines the full range of dispersion of water waves, 
rather than a second order approximation, and the nonlinearity of the shallow 
water equations, and hence it may offer an improvement over the KdV equation for 
short and intermediately long waves. As a matter of fact, numerical experiments 
and 


BKN13 


MKD14 indicate that Whitham’s model approximates waves in 


water on par with or better than the KdV or other shallow water equations do 
in some respects outside the long wavelength regime. (The KdV equation, on the 
contrary, seems a better approximation in the long wavelength regime.) Moreover 
Whitham conjectured that (1.3)-(1.4), where m = cww-: would explain breaking 
and peaking. Wave breaking — bounded solutions with unbounded derivatives — 

|Hurl5| . 


in the Whitham equation for water waves was recently settled in 


___ (See 

|NS94| and [CE98 for related results, and the discussion in [Hurl5|.) Sharp crests 


in its periodic traveling wave of maximum amplitude was numerically supported 
and an analytical result was announced in Ehrl5|. (A complete proof of 


EK13 


global bifurcation will be reported elsewhere.) In the last few years, the Whitham 
equation for water waves received renewed attention from a number of other vantage 
points; see 


EK09 , EGW12 , HJ15 , SKCK14 , EEP15 , among others. 


Benjamin and Eeir in BE67 and, independently, Whitham in Whi67 , on the 


other hand, formally argued that a 27r/A:-periodic traveling wave in water under the 
influence of gravity would be unstable, leading to sidebands growth, namely the 
Benjamin-Feir or modulational instability, provided that kd > 1.363.... Benny 
and Newell in |BN67 


and Zakharov in Zak68 , and later Hasimoto and Ono in 


H072 formally derived a nonlinear Schrodinger equation, which describes the long 


time behavior of the envelope of Stokes waves, and deduced the Benjamin-Feir 
instability. Kawahara in Kaw75 , Djordjevic and Redekopp in [DR77 and later 


Hogan in Hog86 extended the results to capillary-gravity waves. 


tA relative error between cyyw and the phase speed for the KdV equation of, say, 10% is made 
for kd > 1.242 .... 
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The Benjamin-Feir instability is another high frequency phenomenon in water 
waves, which does not manifest in the KdV and other shallow water equations]^ 
Including all the dispersion present in water waves, notwithstanding, the Whitham 
equation may predict short waves’ instabilities. Recently in [HJ15 , as a matter of 
fact, analytically proved was that a small-amplitude, 27r/fc-periodic traveling wave 
of (1.3)-(1.4), where m = cww, be unstable to long wavelength perturbations if 
kd > 1.145..., and spectrally stable to square integrable perturbations otherwise. 
Numerical experiments in SKCK14 bear it out. Note that the critical wave number 


compares reasonably well with that in Whi67 


Here we take matters further and examine the effects of surface tension and 
constant vorticity upon the modulational stability and instability in the Whitham 
equation for water waves. In the presence of effects of surface tension, we extend 
the results in 


HJ15 to (1.3)-(1.4), where 


(1.5) 


{k)=±x{g + Tk^) 


tanh(fcd) 


is the phase velocity for capillary-gravity waves in the irrotational setting; T ^ 0 
is the coefficient of surface tension. At the air-sea interface, q « 9.81m • and 


T « 7.3 X 10 '^N ■ m and thus the effects of surface tension become negligible if 


wavelengths are several times greater than 2tt^T jg ; 
(|1.5|) reduces to (11.21). 


1.7cm. In the case of T = 0, 


We show that the effects of surface tension considerably alters the modulational 


stability and instability of small-amplitude, periodic traveling waves of (1.3)-(1.4) 


and (1.5). In the kd versus k^jT/g plane, specifically, we find several regions of 
modulational stability and instability, whose boundaries are associated with the 
extremum of the group velocity, the resonance between the fundamental mode 
and the second harmonic, the resonance between the limiting long wave and a 
short wave, and a resonance between dispersion and the nonlinearity, which in the 
absence of the effects of surface tension takes place at kd = 1.145 .... See Figure 

and 


Our results by and large reproduce those in Kat76 


DR77 , for instance, via 


formal multiple-scale expansions to the capillary-gravity wave problem. In contrast, 
here we make a rigorous calculation of long wavelengths perturbations to a formal 
approximate model of the problem. These two approaches compare (fortuitously) 
reasonably well. 

When the coefficient of surface tension is large, or equivalently, the fluid depth is 
small so that T/gd^ >1/3, we show that a small-amplitude, periodic traveling wave 
of (1.3)-(1.4) and (1.5) is unstable to long wavelength perturbations, provided that 


the wave number is greater than a critical value, and stable otherwise, similarly 
to the gravity wave setting; see HJ15 . In the case of 0 < T/gd? < 1/3, on the 


contrary, we find three intervals of modulationally stable wave numbers and three 
intervals of unstable wave numbers. 

In the capillary wave setting, i.e. 5 = 0, furthermore, we deduce that a small- 


amplitude, periodic traveling wave of (1.3)-(1.4), (1.5) is modulationally unstable 


regardlessly of the wave number, which is consistent with the results in Kaw75 


* Recently in [HP 15] , a small-amplitude, 27r/A;-periodic traveling wave of the Benjamin-Bona- 
Mahony equation was shown to be modulationally unstable k > y/S, and stable otherwise. But 
the instability mechanism is different from that in the Whitham equation for water waves or the 
water wave problem. 
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and [DR77 , for instance. In the deep water limit as kd —>■ oo, we deduce that 
it is modulationally stable if 1 < Tk'^/g < 1/2, and unstable otherwise, 

which is consistent with the results in Lig65| and [Sini69| , for instance, via formal 
multiple-scale expansions to the physical problem. 


In the presence of the effects of constant vorticity, similarly, we extend the results 


in HJ15 to (1.3)-(1.4), where 


( 1 . 6 ) 


m(k) = 


7 tanh(fc(i) gta,nh{kd) 7 ^tanh^(fc(i) 


2k 


± 


Akd 


is the phase velocity of Stokes waves with the vorticity 7 S K a constant; se^ 


CS04 


for instance. For general non-constant vorticities, the dispersion relation may not 
be written in closed form; see 


HL08 and Karl2 , for instance. 


The zero vorticity setting serves as a reasonable approximation in some circum¬ 
stances and it is easier to handle, analytically and numerically, than the rotational 
counterpart. But the inclusion of vorticities is likely to be far more applicable. As 
a matter of fact, waves are most commonly generated by wind and in any region 
where the wind is blowing there is a surface drift of the water. Moreover currents 
cause shear at the bed of the sea or of a river. Constant vorticity is representative 
when waves are short compared to the length scale of the vorticity distribution; 
lTdSP 88 | 


see 


for instance. 


Colin, Dias and Ghidaglia in CDG95 formally derived an amplitude equation 


for a modulated Stokes wave with vorticity, not necessarily constant, which is the 
same as that in the irrotational setting, whereby they argued that vorticity does not 
qualitatively alter the modulational stability and instability of Stokes waves. (The 
velocity field depends upon the vorticity, though.) But vorticity does influence upon 
the stability of water waves. As a matter of fact, some small-amplitude, periodic 
traveling waves over a shear flow were shown in HL08 to be linearly unstable to 


the same period perturbations, regardlessly of the vorticity strength. Incidentally 
constant vorticity does not seem to invite instability of the kind. 

For each 7 G M we show that a small-amplitude, periodic traveling wave of 


(1.3)-(1.4) and (1.6) is unstable to long wavelength perturbations, provided that 


the wave number is greater than a critical value, and stable otherwise, similar to 
the zero vorticity setting; see HJ15 . But the critical wave number depends upon 


the sign and the strength of the constant vorticity. Specifically, the critical wave 
number increases unboundedly as 7 —>■ —00 in the case of the -I- sign in (1.6) and 
as 7 —>■ - 1-00 in the case of the — sign; see Figur e [4} For each fc > 0, therefore, a 
small-amplitude, 27r/fc-periodic traveling wave of (1.3)-(1.4), (1.6) is modulationally 
stable, for constant vorticities sufficiently large in strength, whose sign depends 


upon the sign in (1.6). Our result helps to shed light into the stabilization effects 


of constant vorticities upon Stokes waves, which recent numerical experiments in 
0SV15 , for instance, suggest. 


^This seminal work quickly led to a flurry of research activities about traveling water waves 
with vorticity, which we do not do justice here. 
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Notation. Let range p G [l,oo] denote the space of 27r-periodic, mea¬ 

surable, real or complex valued functions over K such that 


= I" l/r < +00 if p 


< 00 


and essentially bounded if p = oo. Let denote the space of L|^-functions such 
that f e and 

We express / € d'ln as a Fourier series: 


f{z) 




where /„ = — 


''f{z)dz. 


If / G L 27 I-J -P > 1) then the Fourier series converges to / pointwise almost every- 

1 /2 

where. We define the half-integer Sobolev space H 2 ^ via the norm 


1/2 — /^ + ^ \n\\fn 


n^lj 


2. Small amplitude capillary-gravity waves 

We begin by discussing how one constructs periodic traveling waves of the 
capillarity-modified Whitham equation^ (1.3)-(1.4) and, abusing notation, 

tanhfc 


( 2 . 1 ) 


m{k) = ^gd(l -f 


Moreover we record their small amplitude asymptotics. Note that kd i-G- k turns 
(1.5) to (2.1). One may assume, without loss of generality, the -I- sign in (1.5). As 
a matter of fact, x i— >■ —x and u i— > —u turn one equation to the other. 


In the long wave limit as /c —> 0, one may expand the symbol in (2.1) and write 
that 

m(fc) = v^(l-i(i-^)fc2)+0(fc4). 

Unless T / gdf = 1/3, therefore, the KdV equation with surface tension 


( 2 . 2 ) 


ut + y/gd(l + + = 0 


T 

gd^, 


approximates up to second order the dispersion of the capillarity-modified Whitham 
equation, and hence the water wave problem with surface tension. In the short wave 
limit fc —)■ 00 , on the other hand, m{k) ^ ^JT\k\, and (1.3)-(1.4) and (2.1) become, 
after normalization of parameters, a perturbation of 


Ut 


uux = where A/(fc) = |fc|/(fc), 


in the family of the KdV equations with fractional dispersion, which incidentally is 
L^-critical. 


A traveling wave solution of (1.31 takes the form u(a;, t) = u{x — ct), where c > 0 
and u satisfies by quadrature that 

(2.3) Mu — cu-\- u'^ = (m(0) — c)^b 
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for some b gR. In other words, it propagates at a constant speed c without change 
of shape. We seek solutions of (2.3) of the form 

u{x) = w(z), z = kx, 

where fc > 0 is the wave number and w is 27r-periodic, satisfying that 
(2.4) MkW — cw + w'^ = (to(0) — cYb. 

Here and elsewhere, 

= m(kn)e^'^^ for nGl 

and it is extended by linearity and continuity. 

For each T > 0, note that m is smooth, even, strictly and unboundedly increasing 
away from an interval of low wave numbers. For each T > 0 and k > 0, therefore, 
Aik maps even functions to even functions. Moreover 


Mk ■■ H; 


2-r 


h: 


s+l/2 


27r 

1 


for all fc > 0 for all s ^ 0 


is bounded. Therefore, if ic G satisfies (2.4) for some c > 0, A: > 0 and b gR 
then w G H^. Indeed {Aik — c)w = — + (m(0) — c^b G by the Sobolev 

3 /2 

inequality, and hence w G H2^ . The claim then follows from a bootstrap argument. 
Thanks to Galilean invariance of (2.4) under 

(2.5) w{z)^w{z)+v and ci-Ac—2u, (to( 0) — c)^6 i-A (m(0) — c)^6 —cu+ 


for any u G K, furthermore, we may take, without loss of generality, 6 = 0 in (2.4). 
In the following section, however, we shall make use of variations of the carrier wave 
in the 5-variable. 


For each T > 0 and k > 0 fixed, clearly, the zero function is a solution of (2.4), 
where 6 = 0, for all c > 0. We shall determine at which values of c a branch 
of non-constant even i72.n.-solutions, and hence smooth solutions, emanates from 
the zero solution. It follows from bifurcation theory that a necessary condition is 
that the associated linearized operator admits a nontrivial kernel. While this is 
not in general a sufficient condition, bifurcaiton indeed occurs if the kernel is one 
dimensional. 

A straightforward calculation reveals that the kernel of the linear part Aik — c 
of (2.4) is trivial unless c = m{kn) for some n an integer, in which case 


cos(nz) G ker{Aik — c) 

in the sector of even functions in In the case of n = 1, for instance, bifur¬ 

cation is guaranteed to occur, provided that m{k) ^ m{kN) for all = 2,3,.... 
Otherwise, m{k) = m{kN) for some ^ 2 an integer, and the kernel at c = m{k) 
is two dimensional, generated by cos(z) and cos{Nz). Physically, it means that the 
phase velocities of the fundamental mode and a higher harmonic coincide. In what 
follows, we take n = I and examine when bifurcation occurs at c = m{k). 

When the coefficient of surface tension is large, or equivalently, the fluid depth 
is small so that T/gcP ^ 1/3, it is readily verifed that m(z) > 0 for z > 0 and it 
monotonically increases, unbounded above; see Figure [^a). Therefore m{kN) > 
m{k) for all fc > 0 for all A^ = 2, 3,.... In other words, bifurcation indeed occurs 
from the zero solution at c = m{k). 

In the case of 0 < T/gd? < 1/3, on the contrary, m'(0) = 0, m"(0) < 0 and 
m(z) —>■ oo as z —>■ oo. In particular, the graph of m supports a single local 
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Figure 1. Schematic plots of the symbol m in (2.1) for (a) 
T/gcP > 1/3, (b) 0 < T/g(P < 1/3, and (c) T = 0. 


minimum over (0,oo); see Figure [^b). A straightforward calculation then reveals 
that for each N ^ 2 an integer, one finds a unique > 0 such that m(ktq) = 
m{kNN). Moreover Hn strictly decreases to zero as —>■ oo; see [Mar 13 Lemma 3], 
for instance. In the case of T = 0, furthermore, m{z) > 0 for z > 0, and it 
monotonically decreases to zero as z —>■ oo; see Figure [^c). Therefore m{k) > 
m{kN) for all fc > 0 for all N = 2,3,... and bifurcation occurs at c = m{k); 
see [hjT^ , for instance. 

Let 


( 2 . 6 ) 


E = 


(0,oo) 

(0,oo)\{A:^ :iV = 2,3,...} 


if T/gdP ^1/3, 
if Q <T/g(P <1/2,. 


For each /c S S one may employ bifurcation theory and obtain a one-parameter 
family of smooth and even solutions of (2.4) near the zero solution and c = m{k). 
We summarize the conclusion below. The detail is found in HJ15 , for instance. 


We merely pause to remark that for |6| sufficiently small we appeal to the Galilean 
invariance in (2.5). 

Proposition 2.1 (Existence). Let T > 0. For each k G 'S and for |6| sufficiently 
small, a family of periodic traveling waves of (1.3)- (0, ( |2.1| ) exists and 

u{x, t) = w{a, b){k{x — c{k, a, b)t)) =: w(k, a, b){z) 

for |a| sufficiently small, where w and c depend smoothly upon k, a, b. Moreover w 
is smooth, even and 2'K-periodic in z, and c is even in a. Furthermore, 


(2.7) 


w{k, a, h){z) =WQ{k, b) 
1 


+ 2 ““ 


acos(z) 

' 1 

^m{k) — to ( 0 ) 


cos(2z) 
m{k) — m(2k) / 


0(a(a^ + 6 ^)) 


and 


(2.8) c(k,a,b) =Co(k,b) + (- 

V? 


-k 


m{k) — m(0) 2 m{k) — m(2k) 7 


) -I- 0(a(a^ -I- 6^)) 
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as |a|, l&l —>■ 0, where 

Co { k , b ) = m ( k ) + 2(m(0) — m { k)){b — + 0 { b ^)) 

and 

wo { k , b ) =(to(0) — m { k)){b — b ^ + 0 { b ^)). 

In case k ^ i.e. m{k) = m{kN) for some n ^ TV an integer, the classical theory 
of bifurcation does not apply, but one may nevertheless employ the Lyapunov- 
Schmidt method, for instance, and obtain a branch of small amplitude solutions, 
for instance, for the detail in the capillary-gravity wave problem. 


See Jon89 


In particular, Wilton ripples arise when the phase velocities of the fundamental 
mode and the second harmonic match, i.e. m{k) = m(2k), or equivalently, 

TtVa= . 

3 — tanh^(/cd) 

which in the deep water limit as fed —>■ oo becomes Tk^/g = 1/2. The stability 


theory in the following section, or in HJ15 , does not apply for resonant solutions. 


however, and it is an interesting direction of future investigations. 

3. Modulational Instability Index 

Let T > 0. Let w = w ( k , a, b ) and c = c ( k , a , b ), for fc £ S and |a|, | 6 | sufficiently 
small, form a small-amplitude, 27r/fc-periodic traveling wave of (1.3 )-( |1.4[ ) and (2.1), 
whose existence follows from the previous section. We adapt the developments 
hjTs] (in the gravity wave setting, i.e. T = 0) and study its modulational 


stability and instability. 


Linearizing (1.3) about w in the frame of reference moving at the speed c, we 
arrive at that 

vt + kdz{Mk — c + 2w)v = 0. 

Seeking a solution of the form v{z,t) = A £ C and v £ L^(M), moreover, 

we arrive at the spectral problem 

(3.1) Xv = dz{—M.k + c — 2w)v =: C{k, a, b)v. 

We say that w is spectrally unstable if the L^(IR)-spectrum of C intersects the open, 
right half plane of C and it is (spectrally) stable otherwise. Note that v needs not 
have the same period as w. Since the spectrum of C is symmetric with respect to 
reflections about the real and imaginary axes, w is spectrally unstable if and only 
if the L^(IR)-spectrum of C is not contained in the imaginary axis. 


It follows from Floquet theory (see [BHJ16 , for instance, and references therein) 
that nontrivial solutions of (3.1) cannot be integrable over K. They are at best 
bounded over K, and the L^(IR)-spectrum of C is purely essential. In the case of, 
e.g., the KdV equation, the essential spectrum of the associated linearized operator 
may be related to eigenvalues of a monodromy map, and the stability may be studied 
with the help of Evans function techniques and other ODE methods. Confronted 
with nonlocal operators, ODE methods may not be readily applicable. It follows 
from Eloquet theory (see BHJ16 , for instance, and references therein), instead. 


that any eigenfunction of (|3.1|) takes the form 


j(z) = . 


where (j) is 27r-periodic and ^ £ [—1/2,1/2). 
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Therefore A G C belongs to the (]R)-spectrum of C if and only if 
\v = Cv and v{z + 27r) = 


or equivalently, 

(3.2) \(t> = e-^^^C{k, a, =: C^{k, a, b)(j> 

for some ^ G [—1/2, 1/2) for some 27r-periodic (j). For each ^ G [—1/2,1/2) the 
spectrum of comprises entirely of discrete eigenvalues of finite multiplicities. 
Furthermore 


speCi 2 (ji)(£(A:,a,&)) = |J spec^i^ a, 5)). 

{G[-l/2.1/2) 


Thereby we parametrize the essential L^(IR)-spectrum of C by the one-parameter 
family of point L|^-spectra of £^’s. 

Note that 5 = 0 corresponds to 27r-periodic perturbations, i.e. perturbations with 
the same period as the carrier wave, and | 5 | <SC 1 corresponds to long wavelength 
perturbations of the carrier wave. Slow modulations of the carrier wave are a special 
class of long wavelength perturbations in which the effects of the perturbation 
are to slowly vary, namely modulate, the wave characteristics — k, a, b, in the 
present setting — and the spatial translation mode. Note that variations in these 
parameters provide spectral information about C at the origin in the spectral plane; 
see BHJ16 , for instance. We then say that w is modulationally unstable if the 


spectrum of in the vicinity of the origin is not contained in the imaginary axis 
for some | 5 | small, and it is modulationally stable otherwise. 


The eigenvalue problem (3.2) must in general be investigated numerically, except 
in few special cases, e.g. completely integrable systems. In the case where A is 
near the origin and | 5 | is small, nevertheless, we may take the perturbation theory 
approach in HJ15| , for instance, and address it analytically. Specifically we first 
study the spectrum of the unmodulated operator Cq at the origin; note that zero 
is an eigenvalue of Cq. We then examine how the spectrum near the origin of the 
modulated operator bifurcates from that of Cq for | 5 | small. 

In what follows, thanks to Galilean invariance (see (2.5)), we may take 6 = 0. 


In the case of a = 0, namely the zero solution, a straightforward calculation 
reveals that 


/: 5 (fc, 0 , 0 )e“^ = for all n G Z for all 5 G [-1/2,1/2), 

where = (n -|- ^){m{k) — m{k{n + 5 )))- Notice that 

w±i,o = wo,o = 0. 

As a matter of fact, zero is an L|^-eigenvalue of £o(fc:0, 0) with algebraic and 
geometric multiplicity three, and cosz, sinz and 1 form a basis of the corresponding 
eigenspace. For jaj small, moreover zero is a generalized L\^-e\geuvd.hie of C^ik, a, 0) 
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with algebraic multiplicity three and geometric multiplicity two, and 
1 


Mz) = 


2 (m( 0 ) — m{k)) 

— 1/2 + acos( 22 :) 


= cos 0 + 


1 


m{k) — m{2k) 


{{dbc){daw) - {dac){dbw)){k,a,0)iz) 
0{c 


la^), 


4 ’ 2 {z) = - {dzw){k, a, 0 )(z) = sinz + a- 

a 


sin( 2 z) 


0{a^), 


m{k) — m{2k) 

(l>3{z) =1 

form a basis of the corresponding generalized eigenspace. One may then proceed 


as in 


HJ15 by computing the 3x3 matrix representing the action of on the 
spectral space spanned by (/'s, j = 1 , 2 ,3, associated with three eigenvalues at zero. 


We summarize the conclusion below. The detail is found in HJ15 , for instance, 


which is rather technical, involving commutators of Mk and z. We merely pause to 
remark that an arbitrary to( 0) may be treated by a scaling argument. In the case 
of T/gcP = 1/3, the quadratic order term in the Taylor expansion of m vanishes, 
and the argument in ||HJ15| becomes inconclusive. 

Theorem 3.1 (Modulational instability index). Let T > 0, T/g(P 7 ^ 1/3 and let 


k A 2TT/k-periodic traveling wave of { 1.3 )-( [T4| , (2.1) with sufficiently small 
amplitude is modulationally unstable if IS.mi{ kd) < 0, where 


(3.3) 


A {zm{z))"{{zm{z)y - m(f))) ^ ,, 

^Ml(Z) = -yv-7—y- /SbF\Z) 


m{z) — m( 2 z) 

and Abf{z) = 2(m(z) — m{2z)) + ((zm(z))' — to(0)). It is modulationally stable if 
AMi{kd) > 0 . 

Spectral instabilities away from the origin may be studied using a Krein signature 


method in DT15 , for instance, although we do not consider them here. 


Theorem |3.1| identifies four mechanisms which cause change in the sign of the 
modulational instability index Am/, and hence change in modulational stability 
and instability: 

(1) (zm(z))" = 0 at some z = kd, i.e. the group velocity Cg{k) = {{kd)m{kd)y 
attains an extremum at some wave number k; 

(zm(z))' = m(0) at some z = kd, i.e. the group velocidy coincides with the 
phase velocity Cp(fc) = m(kd) of the limiting long wave at A: = 0, resulting 
in the resonance between long and short waves; 

to(z) = m{2z) at some z = kd, i.e. the phase velocities of the funda¬ 
mental mode and the second harmonic coincide, resulting in the “second 
harmonic resonance”; note from the previous section that k £ T, excludes 
the possibility; 

Abf{z) = 0 at some z = kd, which in the absence of the influence of surface 
tension occurs at kd = 1.146 .. 


( 2 ) 


(3) 


(4) 


see 


HJ15 


Possibilities (l)-(3) are purely dispersive and they do not depending upon the 
specific nonlinearity of the governing equation. As a matter of fact, they agree with 
those in Kaw75| and ||DR77 , for instance, via formal multiple-scale expansions to 


the capillary-gravity wave problem. Furthermore they appear in the modulational 
instability indices for other, nonlinear dispersive equations; see HP15 , for instance. 
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The possibility (4), on the other hand, is associated with a resonance of dispersion 
and the nonlinearity, and it depends upon the nonlinearity of the equation; see 

and 


HP15 , for examples. It does not agree with those in Kaw75 


DR77 . This is 


not surprising since (1.31 does not include many higher order nonlinearities of the 


physical problem. It is interesting to attempt to utilize a more physically realistic 
nonlinearity, but the task is well beyond the scope and interest of the present work. 

The effects of surface tension considerably alter the modulational stability and 
instability in the capillarity-modified Whitham equation. Possibilities (1) through 
(4) may occur depending upon the coefficient of surface tension, the fluid depth 
and the carrier wave number, and the sign of the modulational instability index is 
determined through a rather complicated balance of dispersion and the nonlinearity. 
Figure [^illustrates regions in the kd vers us ky^T /g plane of modulational stability 
and instability. It agrees with those in [Kaw75| and 
branch of the curve along which = 0. 


DR77 , except the “upper” 


In the gravity wave setting, i.e. T = 0, a straightforward calculation reveals that 
the phase velocity m{z) and the group velocity {zm{z)y monotonically decrease 
for all z > 0 with no critical points, and hence possibilities (1) through (3) do not 
occur. Numerical evaluation, on the other hand, indicates that Abf{z) makes one 
sign change, from positive to negative, at z « 1.145... and hence the possibility 
(4) contributes to change in the modulational stability and instability ; see HJ15 


for the detail. Note in passing that it seems difficult to analytically study the sign 
of Abp, but it is explicit, involving hyperbolic functions, and hence amenable of 
numerical evaluation. 

When the coefficient of surface tension is large, or equivalently, the fluid depth 
is small so that T/gd^ >1/3, similarly, the phase velocity and the group velocity 
monotonically increase over the interval (0, oo) with no critical points, and hence 
possibilities (1) through (3) do not occur. Numerical evaluation indicates that 

p makes one sign change, and hence the possibility (4) contributes to change in 
the modulational stability and instability. As a matter of fact, a small-amplitude. 


periodic traveling wave of (1.3)-(1.4), (2.1) is unstable to long wavelength pertur¬ 


bations, provided that the wave number is greater than a critical value, and stable 
otherwise, qualitatively the same as in the gravity wave setting T = 0. 

In the case of 0 < T/gd^ < 1/3, in stark contrast, the group velocity possesses 
a unique minimum. Moreover {zm{z)y = to(0) = y/gd and m{z) = m(2z) each 
takes a root for some z S (0,oo). Therefore each of possibilities (1) through (4) 
contributes to change in the modulational stability and instability, and we find 
intervals of stable and unstable wave numbers. Note in particular that some short 
waves are modulationally stable. Note however that for each 0 < T/gd?' < 1/3, 
small-amplitude, periodic traveling waves become modulationally unstable as the 
wave number increases to infinity. 

In the capillary wave setting, i.e. 5 = 0, we deduce from Figure that small- 
amplitude, periodic traveling waves of (jil -( |1.4[ ), ( |2.1[ ) are modulationally unstable 
for all wave numbers. The result is in agreement with [Kaw75| but it goes against a 
conjecture by Barakat (see Bar , for instance) that capillary waves in very shallow 
water are unlikely to be unstable. 
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Figure 2. Stability diagram in the kd versus k^JTJg plane 
for small-amplitude, periodic traveling waves of the capillarity- 
modified Whitham equation. For T ^ 0 fixed, to interpret, one 
must envision a line through the origin with slope yjrjgd?. “S” 
and “U” denote regions of modulational stability and instability. 
The vertical dashed line represents the positive root of Abf, in 
the absence of the effects of surface tension, approximately 1.146, 
and the dashed line through the origin corresponds to the critical 
value T/gd^ = 1/3. Other curves in the plane represent roots of 
Ami and are labeled according to their mechanism. 


In the deep water limit as kd —>■ oo, note in Figurej^that the curves along which 
m{z) = m{2z) and {znn{z))" = 0 approach the horizontal asymptotes at 1/^/2 

and 2/-\/3 — 1, respectively, whereas the curve along which {zm{z)Y = m(0) and 
the “upper” branch of the curve along which Abf{z) = 0 grow without bounds. 
Therefore possibilities (1) and (3) contribute to change in the modulational stability 
and instability. Specifically, a small-amplitude capillary-gravity wave in deep water 
is expected modulationally stable in the range 2/-\/3 — 1 < Tk"^/g < 1/2, and 
unstable otherwise. The result agrees with that in [Lig65| , for instance, via formal 
argument for the physical problem. The limiting dispersion symbol ^jT\k\ -\- g/\k\ 

fails. 


is not smooth, however, and the present treatment, or that in 


In neighborhoods of boundaries of stable and unstable regions in Figure the 
present treatment fails and one must either go to higher order or consider a different 
method. The resonance between long and short waves was studied in [DR77|, for 
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instance. The second harmonic resonance has been studied by many researchers, 


albeit formally; see, for instance, the book by Craik Cra85 


4. Effects of constant vorticity 

We extend the results of previous sections to the vorticity-modified Whitham 


equation, (1.3)-(1.4), where, abusing notation. 


(4.1) 


mj^±{k) = jd 


tanh(/c) 


tanh(fc) 



■ 


,2tanh^(fc) 
4/c2 • 


As a matter of fact, kd ^ k turns ( 1 . 6 ) to (4.1). Note that ±'m^ ±{k) > 0 for all 
7 S M for all fc > 0. Note moreover that m^^±{k) > mo,±(fc) for all fc > 0 if 7 > 0 
and m^^±(k) < mo^±{k) for all A: > 0 if 7 < 0. In other words, a positive constant 
vorticity increases the phase velocity and a negative constant vorticity decreases. 


One may assume, without loss of generality, the + sign in (4.1). As a matter of 
fact. 


(4.2) 


y^-{k) = m-~f^+{k) for all 7 € K for all k > 0. 


In what follows, we shall take the + sign in (4.1) and suppress it to simplify the 
exposition, unless specified otherwise. In the long wave limit as A; —0, a KdV 
equation whose coefficients depend upon 7 approximates up to second order the 
dispersion relation of the vorticity-modified Whitham equation, and hence the water 
wave problem with constant vorticity. In the short wave limit as A; 00 , on the 
other hand, (I.3)-(I.4) and (4.1) becomes, after normalization of parameters, a 


perturbation of 


Ut 


+ A + UUx = 0, 


where ^f(k) = |A;|/(A:). 

As a matter of fact, it was argued in [Hurl 2 to have relevances to the water wave 
problem in two dimensions in the infinite depths. In particular, it shares in common 
with the physical problem the dispersion relation and scaling symmetry. 

For each 7 G K, note that m{z) is smooth, even and monotonically decreasing 
to zero as 2 ; —)■ 00 ; see Figure Therefore m{k) > m{kN) for all A; > 0 for all 
N = 2,3,.... We then follow the arguments in Section]^ or in HJlSj , to obtain 
for each A: > 0 a smooth, two-parameter family of small-amplitude, 27r/A:-periodic 


waves w = w{k,a,b) of (1.3)-(1.4) and (4.1) traveling at the speed c = c{k,a,b), 


where |a| and | 6 | are sufficiently small. Moreover w and c satisfy, respectively, (2.7) 
and ( 2 . 8 ). 


Furthermore we follow the arguments in Section]^ or in HJ15 , to derive the 


modulational instability index Ami in (3.3), which determines the modulational 


stability and instability of a small-amplitude, periodic traveling wave of (1.3)-(1.4) 
and ( |4l^ . 

For any 7 G K, a straightforward calculation reveals that the phase velocity 
m{z) and the group velocity [zm{z))' monotonically decrease for all 2 ; > 0 without 


critical points, and hence possibilities (1) through (3) following Theorem 3.1 do 
not occur. Numerical evaluation, on the other hand, indicates that Abf{z) 

say, 


m 

and hence 


(3.3) makes one sign change, from positive to negative, at 2 :^( 7 ) 
the possibility (4) following Theorem 3.1 contributes to change in the modulational 
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Figure 3. Schematic plot of the symbols rn^^± for a positive con¬ 
stant vorticity. 


stability and instability. Therefore a small-amplitude, 27r/fc-periodic traveling wave 
of (fOI-dOl) and ( |4.1[ ) is unstable to long wavelength perturbations if kd > Zc, and 
stable otherwise, qualitatively the same as in the zero vorticity setting 7 = 0 . 
Furthermore, a straightforward calculation reveals that zd^) decreases with 7 and 


lim Zcij) = 00 

7—>•—00 


and 


lim Zc{'y) = 0.957.... 

7—)-00 


In the case of the — sign in (4.1), thanks to ( |4.2| ), we obtain the same result. As 
a matter of fact. 

Am/,+ ( 2 ; 7 ) = Am/,-(- 2 , 7 ) 

for all z > 0 for all 7 € K. But 

Abf,-(^;7) = -^BF,+ {z; - 7 ) 

for all 2 > 0 for all 7 S K, and therefore the critical wave number in the case of the 


— sign in (4.1) increases with 7 from 0.957... to infinity. 

To recapitulate, for each 7 G K, a small-amplitude, periodic traveling wave of 


(1.3)-(1.4) and (4.1) is modulationally unstable, provided that the wave number is 


greater than a critical number, and stable otherwise, similarly to the zero vorticity 
setting. The critical wave number monotonically increases to infinity as 7 —> —00 
in the case of the -I- sign in (4.1) and as 7 —>■ -l-oo in the case of the — sign in (4.1), 


and it is bounded below by 0.957.... Figure illustrates regions in the 7 versus 
kd plane of modulational stability and instability. 

In particular, a small-amplitude, periodic traveling wave of the vorticity-modified 
Whitham equation can be made modulationally stable or unstable, by ensuring 
that I 7 I is sufficiently large. For example, consider some kd > 1.145..., so that the 
corresponding, small-amplitude periodic traveling wave of the Whitham equation 
(in the zero vorticity setting) is modulationally unstable. It can be transitioned into 
the region of modulational stability by taking 7 sufficiently large and negative in the 


case of the -I- sign in (^4.l| and positive in the case of the — sign in (4.1). Similarly, 


kd in the range (0.957..., 1.145...), so that the corresponding small-amplitude, 
periodic traveling wave of the Whitham equation is modulationally stable, can be 
transitioned to the region of modulational instability by taking 7 sufficiently large 


and positive in the case of -I- sign in (4.1) and negative in the case of — sign in 


(4.1). Note that all small-amplitude, 27r/fc-periodic traveling waves of (1.3)-(1.4) 


and (4.1) are modulationally stable if 0 < fed < 0.957..., for all 7 G K regardlessly 


of the sign in (4.1). 
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V 


Figure 4. Stability diagram in the 7 versus kd plane for 
small-amplitude, periodic traveling wave of the vorticity-modified 
Whitham equation. The solid and dashed curves correspond to 
the level curves Abf,+{z) = 0 and Abf-{z) = 0 , respectively. 
“S'+” and “S'-” denote regions of modulational stability with the 
corresponding choice of the sign in (4.1), and “[/+ n [/_” is the 
intersection of regions of instability for two equations. Two level 
sets intersect at 7 = 0 and kd = 1.145 .... 


Recent numerical experiments in OSV15 suggest that a small-amplitude, peri¬ 


odic traveling wave in water with constant vorticity can be made modulationally 
stable by ensuring that the constant vorticity is sufficiently large in strength. Our 
results bear it out, observing modulational stabilization in each branch of the phase 
velocity in the vorticity-modified Whitham equation for 7 sufficiently large and with 
an appropriate sign (depending upon the branch). 

Taking the regions of modulational stability in Figure [^together suggests that if 
one were to consider both branches of the phase velocity simultaneously, say, in a bi¬ 
directional Whitham equation, then one would be likely to observe the modulational 
stabilization effects of constant vorticity, independent of the sign of 7 , as reported 
This is an interesting direction for future investigations. 


OSV15 


In the deep water limit as kd —>■ 00 , Figure indicates that small-amplitude, 
periodic traveling waves of (1.3)-(1.4) and (|4.1[) be modulationally unstable for all 


wave numbers for all vort icitie s. In th e zero vorticity setting, this is in agreement 
with the results in BF67 and Whi67 . The limiting symbol g/\k\ is not smooth, 
however, and the present treatment, or that in HJ15 , fails. 
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